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DEVICE AND METHOD FOR CALIBRATING AND IMPROVING THE 
ACCURACY OF BAROMETRIC ALTIMETERS WITH GPS-DERIVED 

ALTITUDES 

5 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention is directed to enhancements for incorporation into 
an electronic navigation device. More particularly, the present invention is directed to an 
10 enhancement of the calibration and hence the accuracy of barometric altimeter 
measurements with the aid of derived altitudes from a global positioning system. 

Description of the Related Art 

In general, altitude measurements are made using two methods of 
15 measurement. One method utilizes a barometric altimeter. Barometric altimeters are 
devices that sense local atmospheric pressure and use a standard model of the atmosphere 
to convert this pressure measurement into altitude. Altitude measurements are referenced 
to height above mean sea level (MSL). 

It is well known that local atmospheric pressure at a given altitude varies 
20 widely due to the effects of weather, solar heating, and other factors. Thus, in order to 
provide an accurate altitude, barometric altimeters must be calibrated to correct for these 
variations. The Global Positioning System (GPS) is a worldwide navigation system that 
can determine a user's position in horizontal and vertical dimensions. However, GPS 
vertical measurements are currently all referenced to the WGS-84 ellipsoid, a purely 
25 mathematical construct that approximates the shape of the earth. The GPS receiver must 
use a model that relates the height above the ellipsoid to the height above mean sea level. 

Further, it is well known that the vertical measurement of a GPS system is 
inherently less accurate than the horizontal measurements. This is due to the fact that 
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GPS satellites are constrained to be above the horizon for signal reception to occur. This 
geometry is less that optimal for measuring the vertical component of a user's location 
simply because there can not be satellites visible below the user (an optimal configuration 
would have satellites above and below the user). All GPS receivers are able to take into 
5 account satellite geometry (Dilution of Precision) and estimates of other satellite-related 
errors (URA) and provide a statistical estimate of the errors in the horizontal and vertical 
measurements. 

In practice, a barometric altimeter typically provides a more stable 
measurement of altitude than GPS over short time periods. However, over long time 

10 periods, pressure variations can be of such magnitude that the barometric altimeter 
measurement of altitude is less accurate than the GPS measurement. The pressure- 
indicated altitude of an uncalibrated barometric altimeter is typically in error by many 
tens of meters due to normal atmospheric pressure fluctuations, weather fronts and other 
sources. However, this error is of a bias like nature - it is slowly varying with time - 

15 resulting in less accurate barometric altimeter readings over long time periods. 
Accordingly, while an altitude determination derived from barometric pressure may be 
meaningfully accurate in a short time frame, over time, the accuracy of such a 
determination becomes undesirable. Conversely, because GPS derived altitude suffers 
from different complementary errors, over a short time period (typically minutes time 

20 frame), GPS altitude measurements are subject to much larger variations than barometric 
altimeter measurements. 

In an attempt to overcome the foregoing, one proposal combines a GPS 
unit and a barometric pressure sensor in the same housing. However, in that proposal the 
pressure sensor is used to augment GPS derived altitude information. 

25 In particular, McBurney et al., in U.S Pat. No. 6,055,477, disclose a 

method of combination or integration of measurements made using two systems to 
provide better availability or accuracy in altitude measurements by estimating a 
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barometric bias using the difference in altitude obtained from the two sources. The 
McBurney et al. method however, fails to recognize that utilizing the difference between 
a GPS derived altitude and a barometric altimeter altitude, as a term in calibrating the 
barometer altimeter, will common mode out any dynamic changes due to movement of 
5 the barometric altimeter and the GPS in tandem. As a result, the McBurney et al. 
approach undesirably requires that the user not change altitude during calibration periods. 

Additionally, in the stated prior approach, an altimeter may only be 
calibrated using GPS derived altitude information when the user is stationary, often 
referred to as a "Calibration mode". The present invention makes no distinction between 

10 "calibration mode" and "navigation mode", indeed the barometric error is constantly 
being estimated and used to calibrate the system. Furthermore, the present invention 
provides a method to statistically determine the need for calibration which results in both 
the calibration and error estimation numbers being calculated and utilized without any 
user intervention, i.e. the user need not place the device in 'calibration mode 7 to obtain 

15 the required bias parameter for calibration computations. 

There exists a need for a method to take advantage of the long term 
stability of the GPS altitude measurement and the short term stability of the barometric 
altimeter measurement to produce an altitude measurement that is stable and accurate 
over long and short time periods. Additionally, the need exists for a method that uses 

20 both GPS-derived altitude and barometric altimeter altitude to produce an altitude 
measurement that is more stable and accurate than either measurement taken alone. The 
need also exists for an improved method to calibrate a barometric altimeter and to 
compute a barometric altitude correction quantity. Particularly, the need exists for a 
method to have GPS and altimeter outputs to be self calibrating while the user is on the 

25 move. The present invention fills the foregoing identified needs, and other needs, while 
overcoming the drawbacks of the prior art. 
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SUMMARY OF THE INVENTION 

It is an object of the present invention to provide an improved navigation 

unit. 

It is another an object of this invention to compute a barometric altitude 
5 correction quantity such that the use of this quantity during calibration results in 
measurements that, are more stable and accurate than GPS or barometric measurements 
taken independently. 

It is another object of the present invention to improve on constraints that 
are required by prior art methods of various types where there is a user distinction 
10 between calibration and navigational modes. 

It is a further object of the invention to provide unique options for 
statistically determining the need for calibration of an altimeter based on discrepancy 
between GPS altitude measurements and other altitude measurements. 

It is another object of the present invention to provide continuous 
15 calibration of an altimeter while the unit is on the move. 

These and other objects are achieved by a portable unit having an internal 
processor. Connected to the processor are, at a minimum, an input (such as a keypad), a 
display, a memory, a barometric pressure sensor, and a GPS receiver, which also 
connects to an antenna, such as a GPS patch antenna. These components, along with a 
20 power supply (such as batteries) are housed within a housing. As will be understood and 
appreciated, the input and display are accessible at an exterior of the housing, in a 
conventional manner. 

A navigation device incorporating the present invention serves as a GPS 
unit, in that GPS signals from a plurality of satellites may be received by the GPS 
25 receiver, such that the processor calculates position information based upon the received 
signals. The conventional use and operation of GPS units is well known, and need not be 
further described. 



4 



Additionally, the present invention addresses an altimeter. In particular, 
the barometric pressure sensor measures barometric pressure and provides the sensed 
barometric pressure information to the processor. The processor, utilizing stored 
software, then converts the measured pressure into an altitude, which may be displayed or 
5 otherwise communicated to the user. The conversion of barometric pressure to altitude 
may be accomplished in any desired and conventional manner. For example, a lookup 
table may be provided in the memory, where the table contains altitude information 
corresponding to known barometric pressures. Thus, an altitude corresponding with a 
sensed barometric pressure may be retrieved from memory and displayed on the display. 
10 Alternatively and preferably, altitude (or elevation) may be calculated using a known 
equation. 

In particular, the present invention provides a unique navigation device 
and method for a navigation device that combines data from a plurality of sensors and 
position information obtained from a GPS, to provide the user with an accurate 
15 representation of altitude information. Additionally, as stated, the simultaneous access to 
GPS information and altimeter information, as well as calculating the difference between 
that information to obtain an indication of bias, allows for features such as automatic 
calibration and calibration while the user is on the move. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

The objects and features of the invention noted above are explained in 
more detail with reference to the drawings, in which like reference numerals denote like 
elements, and in which: 

Fig. 1 is a front view of an illustrative embodiment of a navigation of the 
25 present invention; 

Fig. 2 is an illustrative block diagram of a navigation device that 
incorporates the present invention; 
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Fig. 3 is a graphical representation of a user's elevation trajectory and 
GPS based and barometric altimeter based elevation readings; 

Fig. 4 is a block diagram illustrating barometric altimeter calibration 
processing according to the present invention; and 
5 Fig. 5 is a block diagram illustrating barometric altimeter calibration 

processing according to the present invention. 



DETAILED DESCRIPTION OF THE INVENTION 

With reference initially to Fig. 1, a navigation device that incorporates the 
10 present invention is denoted generally by reference numeral 10. Navigation device 10 
has a housing 12, a display 14, and an input 16, preferably a keypad input. Other known 
inputs, such as a touch screen, may be utilized additionally or alternatively. The housing 
12 is preferably sized to be portable, although the invention is not limited to portable 
units. 

15 With reference to Fig. 2, navigation device 10 has a processor 18. 

Connected to processor 18 are a memory 20, the display 14, the input 16 and a barometric 
pressure sensor 22. Additionally, a GPS receiver 24 is connected to the processor 18. An 
antenna 26, for receiving GPS signals, is connected to the GPS receiver 24. A power 
source, such as batteries, or a battery pack (not shown), is utilized to supply power to the 

20 various electronic components. Additionally, navigation device 10 may include a port, 
such as serial data port, for connecting the device 10 to a remote processor or personal 
computer for uploading information (such as map information) to the device 10, or for 
downloading information (such as route information) to a remote processor or personal 
computer. Alternatively, the device 10 may include wireless communication capabilities, 

25 such that data is received wirelessly from a remote site. As will be understood and 
appreciated, the various electronic components are housed within the housing 12, such 
that display 14 and keypad input 16 are accessible at an exterior of the housing. 
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With reference to Fig. 3, a graphical representation representing a user's 
elevation trajectory (A), a GPS elevation reading (B), and a barometric altimeter reading 
(C) is illustrated. In particular, an exemplary elevation profile of a user using navigation 
device 10 is represented by line A in Fig. 3 In other words, line A represents the 
5 trajectory of a user using navigation device 10 as, for example, he or she travels over 
terrain. The GPS elevation reading is depicted by signal B. As illustrated, the GPS 
elevation reading B varies over time, as the user moves along the path of use. Similarly, 
the barometric altimeter reading varies over time from the actual elevation trajectory of 
the user, although far less significantly, typically, than the GPS elevation reading. 

10 Additionally, as illustrated, the barometric altimeter reading is typically offset by a bias 
amount, or difference D, from the GPS elevation reading. 

With reference now to Figs. 4 and 5, a method of calibrating the altimeter 
of device 10 is illustrated and described. In accordance with an aspect of the invention, 
when the processor determines that the difference between the altitude based upon a 

15 barometric pressure reading from sensor 30 and GPS derived altitude differs by a selected 
threshold amount, the processor begins computations necessary to calibrate the 
barometric readings. In practice, a barometric altimeter typically provides a more stable 
measurement of altitude than GPS over short time periods (from tens of minutes to 
several hours). However, over long time periods, pressure variations can be of such 

20 magnitude that the barometric altimeter measurement of altitude is less accurate that the 
GPS measurement. As such, the processor must determine the appropriate altitude, 
utilizing a combination of these measurements. In particular, as indicated at step 30, the 
device is started up and initialized. Processing advances to step 32, at which the 
processor 18 measures the difference between an elevation reading provided from the 

25 barometric altimeter and an elevation reading provided by the GPS unit. Processing 
advances to step 34, at which the processor 18 computes the average between that 
difference. As will be understood and appreciated, upon the initial measurement at step 
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32, the average difference will simply be equal to the difference. As will be further 
understood, and appreciated, GPS based and barometric altimeter based measurements 
are taken continuously, or periodically, and on additional passes through this processing 
loop, the additional information at step 34 is averaged recursively, although other 
5 averaging techniques may be employed. 

As processing advances to step 36, processor 18 computes the uncertainty 
of the average difference determined at step 34. A decision statistic is employed to make 
the decision to use the estimated barometer altimeter difference to calibrate the baro- 
altimeter reading. 

10 At step 38, the processor 18 determines an average barometer drift for an 

elapsed time associated with the computed average difference determined at step 34. In 

particular, memory 20 preferably has, in table form, the average drift of the barometer 

over time. Following step 38, the processing advances to step 40. At step 40, processor 
18 determines whether the uncertainty of the computed average difference, cr A//w (^), as 

15 determined at step 36, is significantly less than the uncertainty due to baro drift, a Baro> as 

obtained at step 38. In particular, it is determined whether the uncertainty of the 

computed average difference is less than the determined average barometer drift by a 
selected threshold. When the uncertainty of the computed average difference a ^^{k) is 

not less than a Baro by at least the threshold amount, processing returns to step 32, so that 

20 the processor may continue taking difference measurements. When, however, it is 
determined at step 40 that the uncertainty of the computed average difference is less than 
a preselected threshold than the obtained average barometer drift, processing advances to 
step 42. 

At step 42, processor 18 determines whether the elapsed time, that being 
25 the time associated with the measurements taken thus far, is greater than a scaled 
correlation time of GPS vertical errors. When it is determined at step 42 that the elapsed 
time is not greater than a scaled correlation time of GPS vertical errors, processing 
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returns to step 32. When, however, it is determined that the elapsed time is greater than a 

scaled correlation time of GPS vertical errors, processing advances to step 44. At step 

44, the processor 1 8 calibrates the barometric altimeter, pursuant to the processing flow 

of Fig. 5 as discussed below. Following the calibration at step 44, processing advances to 

5 step 46, where the error statistics of the average difference are reinitiated based upon 

calibrated barometric altimeter and processing then returns to step 32. 

With reference particularly to Fig. 5, calibration of the barometric 

altimeter is illustrated and described. At step 48 of Fig. 5, the device 10 of the present 

invention computes a preliminary calibrated barometrically determined elevation. 

10 The calibration process begins with obtaining a calibrated barometric 

altitude H Bcai by subtracting the estimated calibration altitude difference from the current 

barometric altitude i.e. H Bcal (t) = H B (t)- AH caS to remove the bias and thus approximate 

the true altitude H T (t) . This is the approach taken in prior art which is not an optimal 

approach. The present invention goes further than prior art, by using the calibrated 
15 barometric altitude H Bcai to compute a base calibration pressure P Bca! , which is then 

used to compute local altitude H B . 

Processing then advances to step 50, where the processor 18 computes a 
calibrated base pressure value P B . The calibrated base pressure P B is determined by 
solving the following equation for P B , identified as P B ca} . 
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At step 52, processor 1 8 computes a calibrated barometric elevation for 
use in subsequent measurements. In particular, a calibrated barometric elevation, 
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utilizing the computed calibrated base pressure value, P BcaI , is determined according to 
the following equation: 
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H — 
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5 The calibrated barometric elevation is then displayed on display 14 of 

device 10, and used in further processing. 

Accordingly, the present invention employs a method of estimating the 
barometric bias using the difference term Atf^ to common mode out any dynamical 
changes due to movement of the baro-altimeter and the GPS in tandem which is unique 
10 with respect to known methods. In other words, the method of the present invention 

accounts for the fact that changes in altitude by a user are reflected in both the barometric 
altitude reading and the GPS altitude reading, thereby allowing calibration to take place 
while the baro-altimeter and GPS are in motion. A user is not constrained to be 
motionless during the "calibration mode'*. Furthermore, this method allows the 
1 5 barometric error to be continuously estimated and used to calibrate the system when such 
a need is determined by the previously discussed calibration decision model. 

The best known mode for carrying out the present invention utilizes 
models as described below. 
Barometric model 
20 The barometric altitude is modeled by 

H B (0 = H T (t) + B B {t) + Q B (t) (1) 
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where H B (t) is the barometric pressure indicated altitude, H T (t) is the true altitude in 
MSL, B B {t) is a slowly varying bias-like term, and Q B {i) is a zero mean Gaussian noise 
term of variance o 2 Q . Equation (1) shows that indicated barometric altitude is the sum of 

the true altitude, plus a bias-like term that is due to the pressure variation of local 
5 pressure from the standard atmospheric model, and a noise term that is due to noise of the 
sensor, A/D, quantization, and other sources. 

In order to calibrate the barometric altimeter, the bias term B B {t) must be 

determined. 

GPS Model 

10 GPS altitude is modeled by 

H G (t) = H T Q) + B G (t) + C 6 (0 (2) 
where H G (t) is the GPS altitude (in MSL), B G (i) is a slowly varying bias term due to 
ionospheric errors, ephemeris errors, satellite clock errors, and other factors, and C G (/)is 
a zero mean correlated noise term of a much shorter time constant than either B G (t) or 

15 B B {t) . The variance of the C G {t) process is Oy GPS and is an estimate of the errors 

associated with the vertical channel. When Selective Availability was in operation, C G (t) 
was the largest contributor to GPS altitude error (by far). Also, the B G (t) term is typically 
much smaller magnitude than the B B {t) term. 

Determining Error in Barometric Altitude 

20 One approach to calibrating the baro-altimeter using GPS is to simply 

perform a difference of equations (1) and (2) at a particular point in time where certain 
statistical rules (to be discussed later) are met. 

11 



H B (0 - H G (t) = AH(t) = 5, (/) + Q B (0 - 2? c (0 - Q (/) (3) 
Taking the expected value of AH(t) yields 

&H{t) = B B {t)-B c {t) (4) 
since the expected value of terms Q B {t) and C G {t) is zero (they are zero mean noise 
5 processes). 

One can then calibrate the baro-altimeter using AH (/) . Calibration is 
discussed in more detail later. 

Because Q B {i) and C G (t) are zero mean random processes, one can 

reduce the error involved estimating^ (/) by averaging AH (/) over many samples. Note 
10 that B G (t) is ignored since it is typically small. When estimating a random bias in the 
presence of additive noise, the variance of the estimate is reduced by the number of 
samples used to form the estimate only if the additive noise is uncorrected. Q B {t) is 
indeed uncorrected Gaussian noise. However C G {i) is correlated noise with a 
correlation time of x c . Therefore the estimation error is treated differently. 
1 5 One can recursively average AH (t) over many samples according to eqn. 

(5). 

Aff «. (*) = ^Atf^(A-l) + i AH(t n ) (5) 

The uncertainty of this average estimate is the root sum square of 
uncertainty reduction in a Q and the uncertainty reduction to o V GPS eqn. (6). 
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VJGPS 



k*At 



(6) 



where k is the number of samples in the average, and At is the interval between samples. 

Note that the contribution of o V CPS is reduced according to the correlation 

time of this process. What this means is that essentially one correlation time must elapse 
before samples of the C a (t) process are sufficiently decorrelated to contribute an 
uncertainty reduction equivalent to an independent sample. 

It is noted that a y GPS is a dynamically changing function, whereas o Q is 

a quantity that is chosen a-priori. To accommodate these dynamics, g v,gps xs a ' so 
recursively computed over the estimation interval. Again, the same weighting function as 
used in eqn. (5) is used here 




(7) 



Then cr K (&) is substituted into equation (6). 
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Equation (8) is the final form of the uncertainty estimate for AH I 
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Decision Statistics For Calibrating Baro-AItimeter 

The primary decision statistic to use the baro-altimeter difference estimate 
AH^ik) to calibrate the baro-altimeter is when 

where p is any non-negative constant and a Baro (t n - t caJ ) is an estimate of the uncertainty 
to the baro-altimeter. o Baro is a function of the time that has elapsed since the last 
calibration and the uncertainty of the calibration. Furthermore, k is constrained so that 

k*At>a*x c (10) 
which constrains the averaging period to be some multiple a of the correlation time of 

c c (<). 

Also, B B (t) does vary slowly with time. If the averaging period exceeds 
another time threshold, , the assumption that B B {t) is constant does not hold, and the 
averaging process is re-initialized. 

Once the calibration decision statistics have been met, AH^ik) is set 

equal to AH cal9 and all recursive estimation algorithms are re-initialized, in particular 
a Bar0 is set equal toa^^. 

Approach to Calibrating Barometric Altimeter 

The simplest approach to calibration of the baro-altimeter is to simply 
subtract AH caI from the current barometric altitude 

H B , cal (0 = H B (t)-AH CQ! (11) 
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This removes the bias B B (t) and thus H Bcal (t) approximates H T (t) 9 the 

true altitude. 

Indeed, this is the approach taken in McBurney, et al. However, this is not 
the optimal approach. A fundamentally different approach to baro calibration is used in 
5 this invention. 

Standard Atmosphere Model Relating Pressure to Altitude 
For altitude below 1 1,000 meters, the following equation is used to 
compute altitude from pressure. 
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(12) 



10 where the following quantities are define in the 1993 ICAO Standard Atmosphere Model. 
T Q = Standard Temperature at Sea Level 

L = Lapse rate 
R - Gas Constant 
g = Acceleration of Gravity 
15 = Local Pressure (measured by barometer) 

P B = Base pressure (in this case pressure at Sea Level) 
H B = Local pressure altitude 

From eqn. (12), it is shown that the model that relates pressure to altitude 
is an exponential model, not a linear model. Thus, if one determines that AH cal is 30 
20 meters at a nominal altitude of 500 meters, it does not hold that the proper calibration 
factor will still be 30 meters at a nominal altitude of 5000 meters. The reason is that 30 
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meters of elevation difference at 500 meters nominal altitude is a far greater pressure 
differential than 30 meters of elevation difference at 5000 meters nominal altitude. The 
calibration method employed in the invention accounts for this discrepancy. 

Barometric Altimeter Calibration Technique 
In this invention during calibration P L is treated as a constant and P B is 
allowed to vary. The newly computed P B caJ (see eqn. 13) is used in subsequent altitude 

computations in equation (12). 

P B*ot is computed as shown in equation (13). H B cal is the calibrated 

barometric altitude estimated using GPS. 

P B ,,= ^—7 < 13 > 

To 



Summary of Barometric Altimeter Calibration Process 

1 . Evaluate equations (9) and (10) to determine if calibration decision 
statistics are met. 

2. If so, compute H Bcal according to eqn. (1 1). 

3. Compute P Bcal according to eqn. (13). 

4. Begin computing H B according to (12) using P B ca/ . 

5. Set<7 Boro loa^. 

6. Resume estimation of 6H me and Gy^ik) using eqns. (5-7). 
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From the foregoing it will be seen that this invention is one well adapted 
to attain all ends and objects herein above set forth together with the other advantages 
which are obvious and which are inherent to the structure. 

It will be understood that certain features and sub-combinations are of 
5 utility and may be employed without reference to other features and sub-combinations. 
This is contemplated by and is within the scope of the claims. 

Since many possible embodiments may be made of the invention without 
departing from the scope thereof, it is to be understood that all matter herein set forth or 
shown in the accompanying drawings is to be interpreted as illustrative, and not in a 
10 limiting sense. 



